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Abstract: 

We have investigated the modification of physical processes due 
to an intense magnetic field to uncover implications for X-ray spectra 
from binary pulsars. We find that for large enough fields, the quantum 
nature of electron orbits must be taken into account, and for densities 
typical of binary X-ray pulsars, we find that electrons spend most 
of their time in the ground state. For this case, we have estimated 
the production of cyclotron line photons following collisional excitation 
and considered the effects of the transfer of line photons through the 
hot plasma. A line feature as narrow as reported (Trtlmper et al., 1977 
and 1978) appears unlikely. However, more recent observations allow a 
much broader line, and in addition, a consistent absorption feature due 

to knock-on electrons may explain the spectral steepening at ~ 20 keV. 
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I. Introduction 

Pulsars with periods in the range from tens of milliseconds to 
seconds are believed to be spinning neutrino stars with strong magnetic 
fields rigidly anchored to their surfaces. The pulse period corres- 
ponds to the rotation period, and in the case of isolated radio pulsars, 
the rate of change of the period is consistent with the radiative energy 
lost by a spinning magnetic dipole. It. this context, typically, 
the strength of the magnetic field is estimated to be ~ 10 G. 

Fields of this strength alter considerably the interactions of electrons 
with both radiation and ions, primarily by restricting the allowed 
electron recoil transverse to the field. Radiative transitions between 
quantum levels give rise to cyclotron line emission and absorption, and 
the Coulomb scattering process becomes highly modified because of the 
restricted momentum transfer to the electron. The line processes have been treated 
in detail by Daugherty and Ventura (1977 and 1978); the Compton scattering 
cross section for low energy photons (uu « eB/(m e c)) was calculated 
by Canuto, Lodenquai and Ruderman (1971). Resonances in the non-relati- 
vistic Coulomb electron-ion collision frequency were discovered by Ventura 
(1973). Current X-ray observations, however, require more refined 
calculations of the collisional cross sections. In this paper, we 
have calculated an effective scattering cross section for line photons 
and the Coulomb cross section for electron- ion collisions, treating 
the electrons relativistically. 

In section II, we discuss the quantization of electron orbits in 
strong fields, and show that under the, conditions expected at the polar 
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caps of X-ray pulsars, electrons will likely spend most of their time 
in the magnetic ground state, corresponding to zero pitch angle. Then, 
using a fully relativistic quantum electrodynamic approach, we calculate 
the electron- ion Coulomb cross section and the absorption cross sections 
and emission rates for cyclotron photons. Recent evidence has been 
presented for a cyclotron emission line feature at about 60 keV in the 
spectrum of Hercules X-l, an X-ray pulsar in a binary system (Trtlmper 
et al., 1977 and 1978, Coe et al., 1977). In order to explore implications 
of this data, in section HI, we have applied the cross sections developed 
above to the situation where matter accretes onto the magnetic polar cap 
of a neutron star. We calculate an upper limit to the multiplicity of 
cyclotron line photons produced by an accreting proton, and we present 
a rate coefficient for line production in a thermal plasma. Section 
IV consists of an application of our results to the data from Hercules 
X-l. In particular, we consider the production rates of the cyclotron 
emission line mentioned above and its transfer through a hot plasma. 

The line width is estimated as a function of the angle between our line 
of sight and the magnetic field lines, and we find the reported best fit 
line width (Trllmper et al., 1978) places stringent requirements on the 
viewing angle. We also present an interpretation of the spectral 
steepening at ~ 20 keV in terms of an absorption feature which is 
consistent with line observations. Finally, in section V, we summarize 


our results and indicate the need for further observations. 




II. The Effects of Quantization in Intense Magnetic Fields 

The Dirac equation for the spinors ^ representing electrons in 
an external static and homogeneous magnetic field is, in units where 
h * c * 1, 

- m e ) i|r(x) = 0, |i * 0 to 3, (1) 

where m e is the electron mass, are the Dirac matrices, and 

n p " 1 lx^ + e A H (x) ; 


A =0 and "K = % H x x 
o 

The eigenvalues of total energy W are given by 

= p^c^ + m e ^c^ (l+~j), j = 0, 1, 2,.... 


( 2 ) 

(3) 


(4) 


where p is the electron momentum parallel to the field, c is the speed 

is 
2„3 


of light, m e is the mass of the electron, and in cgs units, 


B = m e c ~ 


e h 


4.413 x 10 13 Gauss. 


(5) 


The corresponding spinors for a particular positive energy state, 
assuming B lies in the z-direction, are 

*(x) = 2LW (W + m e ) _J5 (Y^ + + 1PZ U s ; q=<U,2. 

where L is the normalization length,^ is the projection of the position 
vector into a plane perpendicular to the field, s is the spin of the 
electron, and 


u %- 
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0 

0 
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(7) 
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The quantum number j labels the angular momentum parallel to the field 
from the contributions of spin and orbit around the guiding center, 
while q gives the distance from the guiding center to the z-axls of the 
chosen coordinate system. The functions f n ^(x A ) are defined most 
conveniently by the ladder operators: 

~h (tt x + i tt y) and 


then 


and 


/2 

x ~ i 1 ± if ri ^ ! 




X±E „ q ^> " 

together with 


nq' k l' “ 1 eB 1 *n,q » 


f no ( *} = ex P ( - 


eB -*2 


W) 


( 8 ) 

(9) 

CIO) 

(ID 


define the functions completely. The tt operators correspond to the 
canonical momentum of the electron and the X operators to the coordinates 
of the guiding center. There are two contributions to the transverse 
energy quantum number j, expressed by setting 

j = n + s + %, n = 0, 1, 2..., s = +*g (12) 

Then n refers to the Landau level, the state of the orbital motion of 
the electron about its guiding center, and s is the contribution of 
spin to the magnetic moment. Note that only s = -1/2 is permitted for 
the ground state. 
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In a magnetic field, radiation can be emitted and absorbed 

accompanying transitions between transverse energy levels. Under normal 

astrophysical conditions, particle energies are much larger than the 

energy level separation, and the particles consequently have large 

parallel momenta and high quantum numbers. This leads to a classical 

approximation, and T-?hile the radiation is still quantized, the large 

parallel momenta Doppler shift the emitted photons into a continuous 

synchrotron spectrum. However, even for mildly relativistic electrons 

12 

(y < 1.1) in fields greater than 10 gauss, the quantum nature of the 
emission and absorption of radiation must be taken into account. Dis- 
crete cyclotron line emission rates and absorption cross sections 
have been calculated for intense fields by Duagherty and Venture 
(1977 and 1978). We have repeated these calculations in the represen- 
tation above, which is somewhat different from theirs. The QED matrix 
element for a transition is given by 

S fi * ie J' d4x V*> »i (x) (13) 

where the subscripts f and i refer to final and initial states and 
(x) is the photon vector potential, given by 

-. , . .2TT. % +i(iut-Tc*x) - for absorption / 1A v 

V x) = W G 5 + for emission (14) 

where u) and k are the photon frequency and wave-vector, e is the polari- 
zation vector, and V is the normalization volume. The bar over ^ 
indicates the Pauli adjoint of the spinor, \|f^ty°. 

We have confirmed the results of Daugherty and Ventura (1977 and 
1978) for the absorption cross sections and emission rates using the 
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ladder operator representation in cylindrical coordinates. For example, 
the cross section for an electron in the ground state to absorb a fun- 
damental cyclotron photon is given here in the center of parallel momentum 
frame, where the sum of the parallel momenta of the photon and electron 
vanishes: 

-k 


U) cos 3 0) 2 

a a a 


a La) - nl e l ( m a +9w 

a (u),9j U) W ! ( W j+m e ) ( W i+m G ) ( e 

^ S f,5s 2meS f q ( ^ 1+m e^ ^ + W i+«« G J cos 3 0 e 0 + (Wj.-fnig) e x ) 

+6 Sf ^ tu 3 [(Wi+m e + to cos 3 0) 2 E 0 3 + (Wj+m e ) a c,os 3 0 e x 3 ^ 


where 0 is the angle of propagation of the photon and to its frequency, 
and 

W L « m e (1 + (16) 

w i = (">e 2 + uj a cos a Q)^ (17) 

| - -fj Sina 9 (18) 

e and t r are the polarization amplitudes in the ordinary (wave’s 

electric field parallel to B) and extraordinary modes, respectively. 

As we shall now discuss, this cross section is of particular importance 

for intense magnetic, fields, since it will determine the opacity in the 

line feature. We note that equation (15) differs from the classical 
result by terms of order B/B q , which implies that these corrections 

are important wherever the optical depth exceeds B q /B. 

In a strong field, mildly relativistic electrons will spend most 

of their time in the ground state, unless the densities are very high. 

To see this, Figure 1 shows the lifetime x^ of the first excited state 

as a function of field strength for electrons with no parallel velocity. 

This was calculated by numerically integrating over all photon transitions 

to the ground state, according to 


T i - <£/ d3k % 

6 


( 19 ) 
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where (£)d\ is the rate per electron of photon emission into d^k 

centered about 1c and is calculated in the appendix. As mentioned 
before, there are two possible spin orientations for the excited states; 
the lifetime of each is shown separately in the figure. For low 
enough densities, collisional excitation times will be much longer than 
radiative lifetimes; we will return to this after calculating the 
collisional cross sections. 

In a Coulomb collision, the momentum transfer prefers to occur 
perpendicular to the relative velocity. For energetic protons with small 
pitch angle colliding with an electron in the ground state of an intense 
field, the transfer of energy to the electronic plasma is inhibited 
(Basko and Sunyaev 1975). However, under energetically favorable 
conditions, the proton can transfer transverse momentum to the electron 
by excitation of the quantum levels. In addition, head-on or knock-on 
collisions can occur, in which the electrons are not excited, but are 
reflected along their field line. We have calculated the cross sections 
for these processes as a function of the electron energy in the proton's 
rest frame by a single photon exchange approximation. The matrix element 
is given by 

S fi = -ie Vf (x) Y ° A 0 (C) (x) * ± (x) (20) 

(c) 

where A q v ' (x) is the Coulomb potential due to the ion charge. After a 
Fourier decomposition of this potential as outlined in the appendix, we 
find the cross section a ^ for excitation of the jth state from the 
ground state in a Coulomb collision with an ion of charge Ze with zero 
pitch angle is given in the ion’s rest frame by (in cgs units) 
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°i (p> ’ ^ ir f |p= rhv. 


+6 o' 


s '*h 


2B 


j m/cVc 3 } Ej A <P±-^ 


[2W(W+ra e c a )+pc(p^-p)cJ 3 

( 21 ) 
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where r e is the classical electron radius, 2.82 x 10 cm, p and p' 
are respectively the initial and final electron momenta parallel to 
the field, E is the electron kinetic energy, and s’” is its final spin. 
There are two possible values of the final momentum, given by p'* + * 

V" p 2 -2m e ‘ s c 2 B/Bq j, and these are summed over in equation (21) . 
Finally the are defined by 

.,(*>- J i | a - <t - X, (t-x) ;l (22) 

J X t 

The first three are given here in terms of the exponential integral: 
e o (x) - eX£ l (x) 

ei (x) = -1 + (1+x) e X E-^x) (23) 

e 2 (x) = 1+x -(2+x)x e x E^(x). 


The j=o cross section refers to knock-on collision where the parallel 

momentum of the electron changes sign (in the ion’s rest frame), but 

no transverse excitation occurs. Higher j values refer to collisions 

in which the electron is excited to the j level. These results are 

12 12 

shown in Figure 2 for 2 magnetic field strengths, 10 G and 6 x 10 G. 
Because of a breakdown in the approximation, the cross sections become 
infinite at the thresholds of excitation. However, the bandwidth of 
the breakdown is small, so we have arbitrarily smoothed these infinities. 
The approximation breaks down because p , the final momentum, becomes 
small, and the particles are near each other for a long enough time 
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to exchange more than one virtual photon* These Infinities are the 
resonances discovered by Ventura (1973) . As in the non-magnetized case, 
bound states are even possible. The cross section for deexcitation 
from level j to the ground state is given by equation (22) with p and 
p' interchanged on the right hand side. 

From Figure 2, we see that above the threshold, the excitation 
to the j = 1 level dominates, so we can find the critical density below 
which the radiative lifetime of the first excited state is shorter than 
its collisonal deexcitation rate. That is, for densities lower than 
this critical value electrons will spend most of their time in the 
g'tound state. The critical density is defined by 



(24) 


where cr-^ is given by equation (21), by equation (1®), and v is the 

electron velocity. For example, comparing Figures 1 and 2 gives values 

97 12 29 -3 

of ~ 10 'cm (E > 15 keV) in a field of HP Gauss, and ~ 10 cm 

1 ? 

(E > 70 keV) in a field of 6 x 10 ' Gauss. The X-ray emitting regions 
even in accreting binary X-ray pulsars are expected to be much less dense 
than these values. We note also that because the collisional deexcita- 
tion cross section is nearby equal to that for excitation, collisional 
excitation nearly always results in radiative deexcitation. 

Under conditions where the electrons occupy the ground state, 
the opacity in the cyclotron line is determined from the cross section 
given in equation (15) . Since this cross section is much larger than 
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that for Compton mattering, the possibility arises of a source region 
optically thick to line photons but not to the surrounding continumm. 
Since Compton scattering is the most important factor in radiative 
transfer outside the line process, a line feature will appear in this 
case. The width of the line can be roughly determined by finding the 
photon frequency and angle ranges where the optical depth exceeds unity. 
The optical depth is given by 

Tg (tu > 0 ) = J*dl n e J dp f e (p) cr a (u)*,e*) (25) 


where d? is the line of sight increment, f (p) is the electron parallel 
momentum distribution function, and cr is given by equation (15). The 
asterisk denotes the center of parallel momentum frame. If we assume 
that the electron spectrum is uniform throughout the source, the 6- 
function in equation (15) allows us to evaluate the momentum integral 
to obtain 


Wi W+ 

T t) ( co , 0) = N e S f e (p.) S(u>*,e*) — 

J3 + e + 0) ITJ. 


(26) 


w |W + cos0-p + j 

where N g is the column density of electrons in the source, p_^_ and W + 
are the momenta and total energy which satisfy the line resonance condition, 
given below, and S(u)*, 0*) is the right hand side of equation (15) 
without the 6-function. and p + are solved for from 


1 

w + -P+COS0 = 2 ^(— m e 2 - m 2 sin s 0) 
which yields: 

P + = f g ['(l-§) cos0 +((l-g) 2 - !?■ 0 h ], 


(27) 


( 28 ) 


where § is defined by equation (18). We note an important 
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feature: there exists an angular dependent upper bound u) Jnax (0) on the 

line energy where the quantity under the square root in equation (28) 
becomes negative. This limit is given by 

on h 

‘"max sin0 * <29) 

To undertake a study of the radiative transfer of the line photons, 

we note that for n « n , an electron which is excited by absorbing 

ex 

a line photon will deexcite radlatively before anything can happen to 
change its momentum. Thus, the process is effectively a scattering of 
the line photon, and so, in a transfer problem, the photon number is 
conserved. The scattering cross section is the product of the absorp- 
tion cross section and the probability of emission into a given final angle. 
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This gives the expression that should be used at the infinite resonance in 
the Compton scattering cross section of Canuto, Lodenquai, and Ruderman 
(1971). 

III. Cyclotron Line Production in Binary X-ray Pulsars 

The X-ray luminosity of pulsars in binary systems with main sequence 
companions presumably results from accretion. Here we consider the case in 
which the pulsar magnetic field forces plasma to accrete to a polar cap 
with nearly zero pitch angle. The X-rays are generated by the interaction 
of the accretion column falling at its free-fall velocity (~ 0.5c) with 
the background atmospheric plasma. In a steady state situation, the 
atmosphere must be hot enough to avoid the two stream instability, or 
have thermal velocities comparable to the infall velocity, implying 
temperatures corresponding to tens of keV. In other words, when the 
atmospheric density is comparable to that in the accretion column, the 
two stream instability would heat the atmosphere to these temperatures 
much faster than it could cool. 

Under these conditions, if plasma instabilities are not operative, 
one expects the infalling ions to lose energy in Coulomb processes. 

However, since momentum transfer prefers to occur perpendicular to the 
relative velocity vector in a Coulomb encounter, we see that energy losses 
for zero pitch angle protons will be reduced (Basko and Sunyaev' __ 1975) 
because transverse momentum cannot be continuously transferred to the 
background, but must be given up in quanta with (Ap^i ~2m e c B/B^, 
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To Investigate this effect quantitatively, we have used the results 
of the previous section for the Coulomb cross sections cxj in an intense 
magnetic field to calculate the energy loss rate per unit matter traversed 
for a proton with zero pitch angle by Coulomb processes*. 

? J *»P«.<P>v rel (B.P)a j (B.P)AE j (K i p) (30) 

p j »o 

where m is the proton mass, v is the proton velocity, p is the electron 
P 

parallel momentum, f fi is the electron distribution function in p, v re ^ 
is the relative velocity between the particles in an encounter, and AE^ 
is the energy lost in the collision. We assume a relativistic one 
dimensional thermal distribution of electrons, so that f G (p) is given 
by 

f e (p) - exp [ - ( p a +,»/ ) Vi' ] / [ 2n. e K 1 ( ~ e -) ] ^l) 

where K^(x) is the modified Bessel function (Abramowitz and Stegun 1965) 

and T is the temperature in units of energy. The results are shown 

3 2 

in Figure 3 for four temperatures and a field strength of 6 x 10 ‘ G. 

For comparison the energy loss rate in a non-magnetized gas of temperature 

35 keV is shown by the dashed line (Book and All 1975) and in a field 
12 

of 10 G with temperature. 15 keV by the dash-dotted line. For T » 0, 
there is a sharp discontinuity at the threshold for excitation, but 
this disappears when the electrons take on a thermal distribution, and 
for the higher field strength at T «■ 35 keV, dE/dx is relatively constant 
in energy. The energy losses may be approximated with the classical formula 
by replacing the plasma frequency with the cyclotron frequency in the Coulomb 
logarithm. For a neutral magnetized gas, Lee (3.979) has shown that classi- 
cally the cyclotron frequency should replace the excitation frequency in 


this term. 
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Figure 4 shows the Coulomb mean ranges for the temperatures 
discussed above and a field strength of 6 x 10 G. As Basko and 
Sunyaev (1975) pointed out, for high enough energies elastic proton- 
proton nuclear scattering has a mean free path smaller than the Coulomb 
range ( x p _ p e ^ as ~ 50 gm/cm ). In an elastic scattering, the proton 
Is deflected from zero pitch angle, and then Its transverse energy 
can be given up rapidly to ambient electrons, via the excitation of 
plasma waves parallel to the field, for example. Tims , the equation 
for the energetic zero pitch angle proton distribution function is a 
steady state atmosphere can be written 


clIL 


,(E) 


a* & «E» 

where E is the proton energy, x the amount of matter traversed (in 


(32) 


gm/cm for example), $ is the differential (in energy) flux of protons, 
and Opp is the cross section for a nuclear collision with an ambient 
proton, taken from the compilation of Barashenkov and Maltsev (1961) . 

The solution to this equation with the boundary condition that injection 
occurs at the top of the atmosphere, (x « 0) with flux cp Q & (E-E q ) is 

4 <e,x) - ' “Pt-J* 0 af W <»> 


E 


To evaluate the line productivity of the protons in a one 
dimensional thermal gas, we have integrated the flux given by equation 
(33) with the cross sections for knock-on and excitation collisions. 
The production rates at a depth x of electrons in quantum level j with 
momentum p in dp is given by 

<lj(Pi*) • « e j'dp'f e (p' )J E °dE $ (E,x) (p' *)6 (p-p + (p' ,E)) 



(34) 
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where Is the electron number density, $(E,X) is the proton flux given 
by equation (33), v re ^ (E,p') is the relative velocity between the 
particles, v(E) is the velocity of the ion, p'* is the incident momentum 
in the ion rest frame, and the momentum in the delta function is obtained 
from the conservation equations: 

W* = W' *, p* = + - 2eBj , (35) 


where the *'s refer to the ion's rest frame. The total column density 
production rate is found by integrating along dz, the path of the infalling- 
ions : 

N e ,j(p) = S dz 1j^P> x ) = “ f* dx qj(p,x) (36) 


The resulting expression for this quantity, after taking the delta 
functions into account, is 


(37) 


K, i (p> " T % f 0 dE i W I V Ip; <4* ) Vrel ^±’ E) 

The excitation of an electron to the jth level was assumed to 
yield j line photons since the electrons tend to deexcite in single 
step transitions (Canuto and Ventura 1977, Daugherty and Ventura 1977). 

A knock-on (j = 0) was assumed to produce as many line photons as was ener- 
getically possible, by Coulomb collisons with ambient ions. It is 
important to note that here we neglect energy losses of the knock-ons 
to continuum photons by resonant absorption (Pravdo et al., 1978) 
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which results in the conversion of these photons to line photons; this 

process is discussed below. The results for the multiplicity of line 

photons produced under these assumptions are shown for the. four 

12 

temperatures with B = 6 x 10 G in Figure 5a, and Figure 5b shows the 
fraction of the proton's energy put into line photons as a function 
of the Injection energy. For Figure 5b, we assume a line photon has 
an average energy equal to the energy of the 1st excited state with 
P = o. 

In order to investigate the production of cyclotron line photons 
by the high energy tail of a one dimensional thermal electron distribution 
in collisions with ambient protons, we have integrated the cross section 
for the excitation of the first excited state with such a distribution 
function. The resulting line production rate per unit gas density 
squared, is given by 

0^(1) = $ dp f e (p) o^Cp) v(p) (38) 

The results for temperatures of 1 to 100 keV are shown in Figure 6, for 
12 12 

field values of 10 and 6 x 10 Gauss. As can be seen, the line 
production coefficient a, rises quite rapidly as kT approaches TieB/ (m c) 
and flattens somewhat at higher temperatures. 

As noted above, these, results apply to line production by 
electron-ion collisons alone. However, a photon with energy below the. 
cyclotron energy can collide with a fast (f3 ~0.6 to 0.8) knock-on 
electron and be. absorbed as a cyclotron photon because of the Doppler 
shift. When the electron deexcites, the photon tends to be emitted 
preferentially into the electron's forward hemisphere, or downward 
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toward the star’s surface. The importance of this effect may be 
twofold: reemission at angles around 90°, since these have little Doppler 

shift, leads to an increase in line photons; and the absorption results 
in a depletion of the continuum below the line energy, which could make 
a line appear more intense relative to continuum radiation. It is 
difficult to be quantitative about this effect because although the 
production rate of knock-ons can now be calculated, their momentum 
spectrum depends also on the energy loss rate to the radiation, the 
spectrum of which is unknown. It is, however, interesting to note that 
if the field strength is as high as implied by the data of Triimper 
et al. (1977 and 1978), ~ 6 x 10 ' Gauss, then the knock-on process 
is the dominant one for producing line photons under the previous assump- 
tions (neglect of radiation) . So if now we. assume, that the knock-ons 
lose all their energy by converting ~ 30 keV photons to ~ 60 keV 
photons instead of losing 60 keV per photon produced, the the multiplicity 
is only increased by at most a factor of 2. The, knock-on production rates 
are strongly dependent on the. field strength because of the threshold 
effect; at lower field values, relatively more direct excitations than 
knock-ons occur. 

IV. Application to Hercules X-l 

The results of the calculations in the preceding section of the 
fractional energy lost by infalling protons to line photons may be 
compared with recent possible observations of a line feature in the 
spectrum of Hercules X-l at around 60 keV (Triimper et al., 1977 and 
1978, Goe et al., 1977), Assuming that the line is directly produced 
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by the infalling ions and the rest of the ion’s energy goes into continuum 

radiation, the calculations leading to Figure 5 will apply if radiative 

transfer in the source or outside does not convert continuum photons into 

line photons or vice versa. There are two observed components of continuum 

radiation, hard (2 to 20 keV) and soft (.1 to 1 keV) X-rays, with roughly 

37 -1 

comparable luminosities of 1 to 2 x 10 ergs s (Giaccpni et al., 1973, 

Schulman et al., 1975, Catura and Acton 1975), while the observed line 

35 -1 

luminosity is ~ 2 x 10 ergs s . The infalling protons are generally 
assumed to power both continuum components, so at most, only 1/2 to 1 
precent of the energy comes out in the line. Since the effective hard 
X-ray temperature is 20 keV (Pravdo et al,, 1977), we see from Figure 
5b that this fraction is consistent with infall proton energies of 140 
MeV to 180 MeV. However, we still must investigate the effect of 
radiative transfer, in the source and outside. 

To ascertain whether radiative transfer is important, we need to 

estimate the electron density in the source. From Figure 4, we see that 

-2 

the range of accreting matter is at least of order 10 gm cm , corresponding 
24 - -2 

to ~ 6 x 10 e cm . The reported width of the line feature is 12 keV 
(FWHM) , and this places an upper limit on the linear dimension of the 
source parallel to the field, since the line energy varies as field strength. 

_3 

The upper limit is obtained from assuming B « R as for a dipole field: 


1 -i E ^f 


i R t ■ 7]tlo4 “ <fota> 


(39) 


where R is the neutron star radius. Thus the lower limit on the average 
electron density in the source is given by 


^ Ne . , „20 -3/R s 

n e > L * 10 Cm ( 10km ) 


-1 


( 40 ) 


20 


From equations (25) and (15), then, we estimate the opacity to be 

n a e s hc AL . ,R v' 1 . ,T 

T B~ (2nm e c 3 T) 3 5 m e c 3 ^ 50cm * hLOknf * ( '35keV ; 

This implies that if the line of sight covers 50cm or more of the 
source, radiative transfer will be important. 

The preceding considerations apply if the infalling protons 
directly excite the line. If, however, the protons simply thermalize 
without yielding appreciable line emision, the problem becomes that of 
calculating a thermal production coefficient. Applying the results of 
such a calculation shown in Fig. 6 to the observation of Trumper et al. 
(1977 and 1978), we obtain the product of the density squared and the 
volume V of the source region 

n^V =- 2xl0"^cm ^ (42) 

12 

for a temperature of 35 keV and field strength of 6 x 10 G. As before, 
we can set an upper limit on the volume, and thus a lower limit on the 
density, by assuming that the transverse area occupied by the accreting 

/ x ? 

matter is at most tt ’ (0.1R)", or 

A £ 3xl0 10 cm 2 (fo^) 2 (43) 

This limit, together with that expressed by (39), implies 

^ 0 1ri 20 -3,R v“3 

n 3x10 cm (" fo^ (44) 

so that transfer is important in this case, also. Therefore, an 
optically thin treatment of this radiation, as in Daugherty and Ventura 
(1977) is not likely to relate directly to observations. 
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We have considered the problem of line transfer in a rudimentary 

way, by calculating the opacity at a given angle and frequency for a 

thertUal gas of electrons. The motivation for this is the following: 

where the opacity is large, photons are trapped in the source and 

must diffuse out, and at the surface, the intensity radiating outward 

approaches approximately its blackbody value. Therefore, where 

t(u>»0) ^ 1, we assume the line rises out of the continuum. The best 

fit value of the width of the line reported by TrUmper et al.(1977 and 
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1978) is £ 12 keV; as mentioned above, the field implied is ~ 6 x 10 G, 
and in addition, they show evidence that this feature is produced in 
a gas with temperature ~ 35 keV. So for these conditions and four 
values of electron column density, we have estimated the width of a 
cyclotron line feature for several viewing angles 

by finding the frequencies where t = 1. Since Doppler shifting by the 
motion of the electrons is responsible for most of the smearing of the 
opacity, we expect the line to be narrowest at 90° to the field lines 
(and electron motion) and to broaden as the angle departs from this. 
Another factor which is important is the maximum frequency allowed, 
given by equation (29), especially near 90°, where it is responsible 
for the high energy drop in line opacity. 

The results for the line width are shown in Figure 7, for viewing 

o 22 23 24 

angles from 10° to 90 , and electron column densities 10 ,10 ,10 , 
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and 10 cm . As can be seen, for an electron temperature of 35 keV, 

this calculation cannot reproduce the width of the observed feature. 

22 -2 

The column density can't be much lower than 10 cm because of the 
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estimates of lower limits on the density expressed in equations (40) and 
(44). Thus to be consistent with the observation, we are forced to assume 
a lower electron temperature. As mentioned before, lower energy measure- 
ments have placed the temperature greater than 20 keV; so we have shown 
the line width results for this value also. Here we find that the 
reported width can be reproduced if the viewing angle is restricted to be 
very close, to 90° and if the column density is the very minimum allowed. 

We note that the same conclusion about the viewing angle, has been 
arrived at classically (Meszaros 1978) . We also point out that the 
line of sight through the 20 keV source must be less than 100 cm. at 
most. Preliminary results from the A-4 experiment on HEAO-1, however, 
have indicated that if the emission above 40 keV is interpreted as a 
line feature, the best fit width is ~ 60 keV (the lower limit is 
~ 18 keV, D. Gruber, private communication). From Figure 7, we see 
that the viewing angle is restricted to values greater than about 40° 
in this case, and the source column density can be increased considerably. 

Interpretation of the > 40 keV X-radiation as cyclotron line emission 
doesn't explain the sudden steepening of the spectrum at around 20 keV 
(Pravdo et al., 1978). Boldt et al. (1976) proposed an isotropic Compton 
scattering as the cause; however, we have attempted to attribute this feature 
to an absorption edge or dip. As has been shown , cyclotron absorption can be 
treated as a resonance in the Compton process. A problem arises if the dip 
is due to cyclotron absorption, however, because of the high field strength 
implied 2 x 10 x G) . If the absorbing electrons are located in such a 
high field, then they must be near the hard X-ray source. Consequently, 
they are bathed in X-rays and should be heated to a temperature approaching 
that of the source region itself. We have shown above that line transfer 
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is equivalent to scattering with a much smaller mean free path than 
Compton scattering. In this case, electrons near the source heated 
to temperatures near the source temperature become part of the very 
thick line scattering region. In other words, they re-emlt the absorbed 
radiation in a manner similar to the source itself, and no absorption 
dip appears. The knock-ons discussed in the preceding section are 
a possible solution to this problem. These electrons stream down the 
field lines toward the stellar surface with $ ^ 0.6 to 0.8, and can 
consequently absorb upward moving photons with frequencies less than 
the fundamental cyclotron frequency because of Doppler blue-shifting. 

Since the knock-ons are mildly relativistic, they preferentially reemlt 
a line photon in their forward direction, i,e. downward toward the 
surface out of our viewing direction. If the column density of knock- 
ons is greater than unity, photons with high enough energies could be 
depleted in the same manner as for ionizing photons in an ionization 
bounded HII region. That is, photons Doppler shifted to cyclotron 
frequencies would be screened out by the knock-ons in the same way that 
hydrogen atoms screen ionizirtg photons. 

The absorption depth of the knock-ons can be estimated by using 

the production rate expressed by equation (39) and by assuming that 

the fast electrons lose energy by absorbing upward photons and re-emitting 

them downward. We estimate the accretion rate by assuming an efficiency 

of 10% for the conversion of kinetic energy to X-radiation, and the 

10 2 

cross-sectional area of 3 x 10 cm for the accretion column. The life- 
time is roughly given by the number of photon collisons required to 
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give up its energy 10 for a 250 keV knock-on) divided by the 
collision frequency. The photon number density is obtained from the 
luminosity L , greater than 20 keV, by 

A 


n 


X 


L x (>20 keV) 
A c 


/ , a 21 -3 

4 x 10 cm , 


(45) 


assuming 1 % of the total luminosity comes out above 20 keV. For infall 
energies of VL50 MeV/nucleon and B * 6 x 10^G, equation (37) yields 

t 

N e o ^ % and thus the absorption depth r is estimated to be; 


T a ~ % T ~ cp Q cr a 


10 

n xPa c 


1 


(46) 


for the assumptions given above. This value implies that knock-ons could 
be responsible for an absorption dip. 

This picture also predicts a certain behavior of the absorption 
feature with pulse phase. In particular, as the line of sight comes 
closer to the polar field lines, the energy at which absorption begins 
should become lower, and the variation in energy should become smoother. 
In recent paper (Pravdo et al., 1978), we assumed an arbitrary high 
energy electron distribution, and directly compared data from the cosmic 
X~ray experiment on 0S0-8 with these predictions. As discussed there, 
the fit was acceptable over the region of the spectral pulse, and the 
best fit angles generally performed ^20° excursion in a region around 
40°. As opposed to the line emission interpretation, however, here 
the spectral pulse (which may only be loosely correlated with the 
intensity pulse) defines a pencil beam 
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(or cone). The best fit field strength was coincidentally consistent 

12 

with the emission measurements, ^6 x 10 G. The cyclotron photon 
energies here are too broad to form a well-defined narrow line feature, 
but can be made consistent with the HEAO-1 A4 data. 

V. Summary and Discussion 

12 

The strong magnetic fields 10 G) believed to exist at the 
surfaces of neutron stars can strongly modify the interactions of 
electrons with ions and radiation. We have considered the effects of 
such modifications on the X-ray spectra of the pulsars in binary systems 
which derive their luminosity by accretion from their companions. The 
effect of the intense field is to quantize the allowed electron momenta 
transverse to the field lines. Radiative transitions between levels 
give rise to cyclotron emission and absorption, and since transverse 
momentum transfer from an ion to an electron must occur in discrete 
quanta, the Coulomb cross section becomes a knock-on cross section and 
a series of excitation cross sections. We have shown that electrons 
in the polar caps should almost always be found in their ground state 
and have then calculated the cross sections for these processes treating 
the electrons relativistically. 

Since Hercules X-l is the best studied of these binary X-ray 
pulsars and shows a possible line feature, we have used the above 
cross sections in an attempt to understand the physical conditions in 
the polar cap of this source. Data from 2 to ^20 keV imply that the 
effective electron temperature in the source is greater than 20 keV. 

To produce the observed line feature with a width (FWHM) less than 12keV,we find 
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that the temperature cannot be higher than 20 keV, the angle between 
our line of sight and the polar cap field lines must lie within 10° of 
90® during the entire high intensity part of the pulse, and the linear 
extent of the source along the line of sight must be less than 100cm. 
Recent measurements, however, relax these requirements considerably, 
allowing the viewing angle to be as low as ^40® and the line of sight 
to be much larger. Also, we have attempted to explain the spectral 
steepening at ^20 keV in terms of an absorption dip. The absorption dip 
can be made consistent with line observations because of the ability 
of the ions to produce a high parallel momentum population of knock-ons, 
which would absorb upward moving continuum photons and reemit them 
preferentially downward. 

The observations required to determine the nature of the Her X-l 
radiation, as well as that from similar sources, must include both 
spectral and temporal resolution. The 0S0-8 observations of Her X-l 
(Pravdo et al. 1977 and 1978) showed dramatic spectral changes below 
30 keV over a portion of the pulse phase, and recent observations from 
the A-2 experiment on HEAO-1 (Pravdo et al. 1979) show evidence for 
pulse phase variations at higher energies. Generally pulse periods 
are ~ a second or- longer, so time resolution on the order of tens of 
milliseconds is adequate. As for spectral resolution, we note the 
disagreement between scintillation spectrometers for the higher energy 
radiation, and point out that solid state devices could be made sensitive 
to energies as low as 40 keV with much improved resolution. Perhaps 
most importantly, an understanding of the variations already observed 
at lower energies with proportional counters may make clear the physical 
processes and conditions in such sources . 

For helpful discussions, I would like to thank Drs. Steven Pravdo, 
Reuven Raraaty, Elihu Boldt, Barham Smith, and Joseph Daugherty. Also, 
Dr. S. Langer is gratefully acknowledged for an important correction. 
The research was supported by NASA grant NGR 21-002-316 • 


REFERENCES 


Abramowitz, B. and Stegun, I. A. 1965, Handbook of Mathematical Functions 
(New York; Dover). 

Barashenkov, V. S. and Maltsev, V. M. 1961, Fortschrltte der Physik, 

9, 549. 

Basko, M. M. and Sunyaev, R. A. 1975, Astron. and Astrophys, 42, 311. 

Bjorken, J. D. and Drell, S. D. 1964, Relativistic Quantum Mechanics 
(New York; McGraw-Hill). 

Boldt, E. A., Holt, S. S., Rothschild, R. E. , and Serlemitsos, P. J. 

1976, Astron. and Astrophys. 50, 161. 

Book, D. L. and All, A. W. 1975, A Collection of Plasma Physics Formulas 
and Data, NRL Memorandum Report 2898. 

Canuto, V., Lodenquai, J. , and Ruderman, M. 1971, Phys. Rev. D, 3, 2303. 

Canuto, V. and Ventura, J. 1977, Fundamentals of Cosmic Physics, 2, 203. 

Catura, R. C. , and Acton, L. W. 1975, Ap. J. Letters, 202, L5. 

Coe, M. J. , Engel, A. R. , Quenby, J. J. , and Dyer, C. S. 1977, Natue, 

268, 508. 

Danese, A. E. 1965, Advanced Calculus, Vol. I (Boston: Allyn and Bacon). 

Daugherty, J. K. , and Ventura, J. 1977, Astron. and Astrophys., 61, 723. 

Daugherty, J. K. , and Ventura, J. 1978, Phys. Rev. D, 18, 1053. 

Giacconi, R. , Gursky, H. , Kellogg, E. , Levinson, R. , Schreier, E., and 
Tananbaum, H. 1973, Ap. J., 184, 227. 

Lee, J. F. 1979, preprint. 

Merzbacher, E. 1970, Quantum Mechanics (New York; Wiley). 

Meszaros, P. 1978, preprint. 

Pravdo, S. H. , Boldt, E. A., Holt, S. S. and Serlemitsos, P. J. 1977, 


Ap. J. Letters, 216, L23. 


Pravdo, S. H., Bussard, R. W. , Becker, R. H., Bo. I cl t, E. A., Holt, S. S., 
Serlemitsos, P. J. and Swank, J. H. , Ap. J., 225, 988. 

Pravdo, S. H., Bussard, R. W. , and White, N. E. 1979, M.N.R.A.S., in press 
Schulman, S., Friedman, H., Fritz, G. , Henry, R, C, , and Yentis, D. J. 
1975, Ap. J. Letters, 119, L101. 

Trumper, J. , Pletsch, W. , Reppln, C., Sacco, B., Kendzlorra, E., and 

Staubert, R. 1977, Texas Symp. on Relativistic Astrophys. (Boston, 1976) 
aim. N.Y. Acad, Scl., In press. 

Trumper, J., Pletsch, W. , Reppln, C. , Voges, W. , Staubert, R. , and 
Kendzlorra, E. 1978, Ap. J. Letters, 219, L105. 

Ventura, J. 1973, Phys. Rev. A, 8, 3021. 


Figure Captions 


Figure 1 - The radiative lifetime of the first excited magnetic state 
as a function of field strength, for both Landau (s«-Lj) and 
spin (s*+^) excitations. 

Figure 2 - The Coulomb cross sections for excitation and knock-on 

12 i o 

collisions in fields of 10 and 6x10 G as a function 
of the electron energy in the proton rest frame or vice 
versa, The electrons are initially in the ground state, 
and the Ions have zero pitch angle; j labels the final 
transverse momentum state. 

Figure 3 - The Coulomb energy losses suffered by zero pitch angle. 

protons per unit matter traversed for several temperatures 
in a field of 6x10^ 2 G as a function of the proton energy. 
Also shown are the losses in a non-magnetlzed gas with 
kT » 35 keV and in a field of 10 12 G with kT » 15 keV. 

Figure 4 - The Coulomb range of zero pitch angle protons, neglecting 
nuclear encounters, in gases at several temperature and 
B » 6xl0 12 G. 

Figure 5 - a) The. multiplicity of line photons produced as a result 
of Coulomb interactions with infalling protons as a 
function of the initial energy of the proton, 
b) The fraction of the proton kinetic energy going into 
line radiation. 

Figure 6 - The. coefficient for thermal excitation of the first excited 
state in fields of 1 and 6xl0^ 2 G as functions of the gas 
temperature. 

Figure 7 - The width of the fundamental cyclotron line determined from 
opacity (self-absorption) for two temperatures and four 
values of column density for the source electrons. 
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APPENDIX 


In this appendix, we present the details of the calculations of the 
cross sections and emission rates presented in this paper. First we eval- 
uate an integral occurring in all the calculations: 





in an operator representation. The argument of the exponential can be 
written using aquation (8), 



( k.‘* + ) = ;K + (X'+^r r)+u-(X + -^ir») (AiO 


where a vector V is decomposed according to 

v ± = s: i (as) 


The commutation relations among the various operators are as follows: 

{ if*, ir*3_ = -eB 

CX-, -y-l_ - ^ (ft*) 

[tr 1 , X 1 ]. = o 

Prom Msrzbacher, (1970, p. 167) we obtain for this situation: 


** |= (a, + A 2 ) = exp(^) , ex^c^) . exp (- i ) 


(As) 


Applying this relation to the exponential in A-l gives: 

^ = e *]P^ u * + X ) ex jpL <■ X + ) exp( — ) exjn(- 

ti 

where §- k t ^/(2aB). Thus, substituting in A- 2 and using the hermi'city of 
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rr&aa a sa 


Ch^ operators, (X")+ « X+; (rT) * * 


rr.‘ , yields 




• ^|o(ife.-A + )xxj.(- ■£*-) , 

expanding the exponential and operating on the functions according to 

equations (9) and (10) gives 


(A 7) 


-p-3!ji y F«r' (-Lk)^ Ck. s- 
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.-J 
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which yields after using the orthogonality property of the f's and 

*£* 4, t ™ 

k~ = k ± e“ “ /, 2 , where r is the axinuth of t in the transverse plana, 
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Substituting j! 


^1 > n i -n 2— anc * M-l - max (o , yields 
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Fron the definition of the Laguerre polynominal: 
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we see that (A-10) can be written 

T _ -S * (lV - ^*t-Hn,~nil-vr\ a -n,) r f n\.n(n,,nQ>. w.i.nCg,,^0\ l£ 
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We can now proceed to the evaluation of the cross sections. 

First we calculate the absorption cross section and emission - 
rate for line photons as prescribed by equations (13) and (14) : 

S f \ =!.«. (uv-) 4 ■$<*> ^ mO , (A-ii) 

for absorption. Substitution from equation (6) for % gives 

S rl = XL l' y. (\J f + C *^ 


which, after using the hermiticity of the Dirac matrices (y^ = y^^ 0 ). 

r* ?. e (k> v «.) W s; w 11 ***# 


This expression can be reduced to 2x2 from 4x4 by using 

u s - I'D , »•. U Uf) 

where I s > » (^®*” S , ) and a are the Pauli matrices. After commuting 
«s»-* . 

the exponential factor to the left, we obtain the matrix operator 


(A-lk) 


between states 


y® [ ( ri+kJ+Me T *• s (» M n , y+m a - 


(^-17) 


( (jr© fv,, + ,r O Ct*TT +1 fl“. £xTr } v f-n iff 


where some components are not shown because they are not needed. Then 
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the matrix element A-15 becomes 


c ic r air 1 '4. f r * \ 

S tL " il . U V W t WilW^'KWi+w.i ] J * *•*>«. V •‘ e 

This equation is valid £or all transitions between levels. We 
simplify by requiring that the initial state be the ground state, by 
letting q^ ■ 0 without loss of generality, and by evaluating in the 
center of parallel momentum frame, where pf - o and Wf - W^ 3 /m y -h2eB' 
Under these conditions, equation A-18 becomes 

c I £ C^l £ c ps') [ { 

+ £*-k*£")] 

+ ^ [-K+^KOf-iceO - I^. 0> <> Ux ) . 

We note that if 9 is the polar angle of k and cp its azimuth, 

e*‘ f U.COS& Tle«) C A-io) 

where e Q and e x are the polarization amplitudes in the ordinary and 
extraordinary modes. Substituting for the I’s from A-12 gives 

tafffe $Wt +o-Vt)S ( ‘h'i p-s-i^ (ifi $}*?[( r . /r-ff 
t'u" it j-uw, w; tVi 

[ (h/£+W s **05efls6£o +C^ + ,r *\!)(c*s8foH€ x 5)fS, #> A w[W; +>rj < . +wcas v ^')^o 

4*t ^ €x J J (a 




Carrying out tha same analysis for the emission case gives (assuming 
pi' Si a 0) 

c v Uff/ 4 £ (V, -ui-Wf) S(h« + k*) -3+tjP 

(cnus ^“a.L [uVW,(,W,*m e )W # (*^«.«)li 1/^T 

^ 4r> v)U»iS , S ft -*-U x )l + & s , ^ U3 [(w + ♦ w e +u>C0S*fi) 

+1 W|+w e ^ f o* ^ , 

The absorption cross section is obtained by squaring the magnitude 
of Sj£ and sternal eg over electron final states. The square of the 
o-functions are handled hauristically in Bjorken and Drell (1964) : 

i«c»r- biv , u (vn r = & ^w) , (a») 

where L and I are the (large) normalization length parallel to the 
field and tha normalization tine. Then the rate of absorption is 

. 7 c ^ ^ h l tvi* ■(*«> 

v-° . ... ‘ 

where V is the electron normalization volume, cr a is the absorption 
cross section, and v re ^ is the relative velocity of the photon and 
electron. Solving for tha absorption cross section yields equation 
(15) in the center of parallel momentum frame. 

The emission rate is shown in here in more detail to show the 


processes involved in arriving' at equation (15). First, the sum over 



The integral over is ignored because of tha 6 function. The the 
emission rate into d-*k cantered on k for electrons deexciting from 


the flr3t excited state can be written as 


<- '"MVi - ^ ^ i I vi\ 
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where Vd\/(2n)^ expresses the proportionality of the rate to the phase 
space volume for the emitted photon. Carrying out the operations yields 
„ ) Jit 5 c" T f F r , x _ , 1* 
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(/t-J.7) 
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When carrying out the Integral over d J k to obtain the inverse of the 
lifetime r^* it must be noted that depends on u> so that the 5 func- 
tion in A-27 should be written 
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where the braces represent the quantity in braces in A-27. This integral 
was evaluated numerically; the results are shown in Figure 1. When 
the radiative lifetime of an electron in the first excited state* is 
short compared to collisional timescales, we can combine the absorption 
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and. emission processes into a scattering cross section, givenby 

<L<rCv,p ) tr v e .** K<*>- • (W r ij) (m* faWi^Q ^ g 
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where everything is as defined for equation (15) in the text, and 
the primes indicate the emitted photon. The T^'s are shown in Figure 1. 
We note the symmetry between final and initial photons. The emitted 
photon has the energy 
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Next we consider the scattering of an electron by an ion in an 
external field, first order approximation, the relativistic. extension 
of the Born approximation. The matrix element involved in the ion's 
rest frame is given by equation (20); we Fourier expand the ion potential 


-r/r 


Ze/r-e sc, where r c , is the screening length: 
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Then equation (20) becomes 
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Substituting for from equation (6) and using he.rmiticity gives 


5 A' J k 4~h, Jet 4 * 
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* ik.x 

Arter moving e to the left, we reduce the inner to the 2x2 
representation as in A-16 through A-18: 

U 3 J [* 9 TV if. (TT + k) + *l t ] »* UX +rn e ) U-s L (A3S) 

= | [ (TT 0 + m e ') X + . T«(ir+^ )“«• t'ir+^)xTT ] | , 

Va note that tt x tt * ■£ 7 x (sA) = -r B. Then substituting this expression 
into A-3i and assuming the initial state is the ground state gives: 
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The elasticity of the interaction is expressed by the fact that 
the Coulomb potential had nc time varying component. This results in 
a factor 5 r-W.,. ) in the above tine integration. In evaluating the inte- 
gral, we put vT,. = v*' c - W to obtain: 
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Substituting for the I function and carrying out the kj r and cp integra- 
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IX 


geven area A(= ttR^^ ) is given by 
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Tna rate is obtained from 
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In the squata, the 5 function is handled as in A-23. The 6-function 


remains, and is used in tha p f integral, after 
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Solving (A-42) for cr^ (noting that v = p/w)givesl 
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The sum can be 


reqritten by using first of ill t 
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